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(54) Method for extracting a three-dimensional model from a sequence of images 

(57) In a computerized method, a three-dimensional 
model is extracted from a sequence of imag3? ihat 
includes a reference image. Each image if. the 
sequence is registered with the reference image to 
determine image features. The image features are used 
to recover structure and motion parameters using geo- 
metric constraints in the form of a wireframe mode. A 
predicted appearance is generated for each image 
using the recovered structure and motion parameters, 
and each predicted appearance is registered with the 
corresponding image. The recovering, generating, and 
registering steps are repeated until the average pixel 
value difference (color or intensity) between the pre- 
dicted appearances and the corresponding imao«s is 
less than a predetermined threshold. 
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Description 

Field Of The Invention 

5 [0001] This invention relates generally to computer vision, and more particularly to recovering structure from motion 
as expressed by optical flow in a sequence of images. 

Background Of The Invention 

10 [0002] In computer vision, one problem is to recover the three-dimensional configuration of an object from a sequence 
of two-dimensional images acquired by a camera. This is :</^y difficult when both the camera parameters and 
point correspondences across the image sequence are unknown. 

[0003] There is a large body of work on the recovery of raw 3-D data from multiple images; they include multibaseline 
stereo, trinocular stereo that combines constant brightness constraint with triiinear tensor, stereo with interpolation, and 
15 shape from rotation. 

[0004] Virtually all stereo approaches assume a fixed disparity throughout once the disparity has been established, 
e.g., through a separate feature tracker or image registration technique. Most techniques assume that the camera 
parameters, intrinsic and extrinsic, are known. For 3-D facial modeling, the following techniques are generally known. 

20 From range data: 

[0005] Range acquisition equipment includes light-stripe rangefinders, and laser rangefinders. Rangefinders, when 
compared to video cameras, are relatively expensive, and considerable post-processing is still required. For example in 
one method, feature-based matching for facial features, such eis the nose, chin, ears, eyes, are applied to dense 3-D 
25 data to initialize an adaptable facial mesh. Subsequently, a dynamic model ol facial tissue controlled by facial muscles 
is generated. In another method, a range image with a corresponding color image of a face is used. The 2-D color 
image is used to locate eyes, eyebrows, and mouth. Edges in color space are determined, and contour smoothing is 
achieved by dilation and shrinking. 

30 From two 2-D images: 

[0006] Two orthogonal views of a face are normally used. The profiles are extracted and analyzed; this is followed by 
facial feature extraction. A 3-D face template is then adjusted by interpolation, based on the extracted information. 

35 From a sequence of temporally related 2-D images: 

[0007] In one approach, 2-D images are used to reconstruct both shape and reflectance properties of surfaces from 
multiple images. The surface shape is initialized by conventional stereo image processing. An objective function uses 
the weighted sum of stereo, shading, and smoothness constraints. The combination of weights depends on local tex- 

40 ture, favoring stereo for high texture with a known light source direction and known camera parameters. 

[0008] A calibrated stereo pair of images has also been used. There, a disparity map is determined, followed by inter- 
polation. In one implementation, three-dimensional deformation is guided by differential features that have high curva- 
ture values, for example, thB nose, and eye orbits. If the motion between images in a sequence is small, then the optical 
flow can be used to move and deform a face model to track facial expressions. Fixed point correspondences are defined 

45 by the optical flow. The deformation of the face model is constrained and specific to faces. Facial anthropometric data 
are used to limit facial model deformations in initialization and during tracking with the camera's focal length approxi- 
mately known. 

[0009] In a different approach, facial features such as the eyes, nose and mouth are tracked using recursive Kalman 
filtering to estimate structure and motion. The filter output is used to deform the shape of the face subject to predefined 
so constraints specified by a linear subspace of eigenvectors. 

SUMMARY OF THE INVENTION 

[001 0] Provided is a computerized method for recovering 3-D models from a sequence of uncalibrated images with 
55 unknown point correspondences. To that end. tracking, structure from motion with geometric constraints, and use of 
deformable 3-D models are integrated in a single framework. The key to making the recovery method is the use of 
appearance-based model matching and refinement 

[001 1] This appearance-based structure from motion approach is especially useful in recovering shapes of objects 



2 



EP 0 907144 A2 



whose general structure is known but which may have little discernable texture in significant parts of their surfaces. 
[0012] The invention, in its broad form, resides in a computerized method for extracting a three-dimensional configu- 
ration of an object from two-dimensional images, as recited in claim 1. 

[001 3] The method can be applied to 3-D face modeling from multiple images to create new 3-D faces for a synthetic 
5 talking head. The talking head includes a collection of 3-D triangular facets, with nodes as vertices. The model can be 
recovered even when the sequence of images are taken with a camera with unknown camera focal length and extrinsic 
parameters, i.e., the camera pose is unknown. 

[001 4J In the general method, a three-dimensional configuration is extracted from a sequence of images including a 
reference image. Each image in the sequence is registered with the reference image to determine image features. The 
io image features are used to recover structure and motion parameters using geometric constraints on a 3D wireframe 
mesh template. 

[0015] A predicted appearance is generated for each image using the recovered structure and motion parameters, 
and each predicted appearance is then registered vAlh the corresponding image. The structure and motion recovery, 
appearance generation, and image registration steps are repeated until a selected termination condition is reached, for 
is example! an average difference between the predicted appearances and the corresponding images is less than a pre- 
determined threshold, or a fixed number of iterations has been performed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 [0016] A more detailed understanding of the invention may be had from the following description of a preferred 
embodiment, given by way of example, and to be understood with reference to the accompanying drawing wherein: 

♦ Figure 1 is a block diagram of a 3-D structure recovery system that uses the invention: 

25 ♦ Figure 2 is a flow diagram of a recovery method according to a preferred embodiment of the invention; 

♦ Figure 3 is a diagram of spline control grid superimposed on pixel locations; and 

♦ Figure 4 is shows a facial wireframe model superimposed on an image of a face. 

30 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 
Introduction 

35 [0017] Described is a computerized method for recovering 3-D models from a sequence of uncalibrated images. with 
unknown point correspondences. In order to perform the recovery, the difficult point correspondence and occlusion 
problems must be overcome. To that end, tracking, structure from motion with geometric constraints, and use of deform- 
able 3-D models are integrated in a single framework. The Key which makes the method work is the use of appearance- 
based model matching and refinement 

40 [0018] This appearance-based structure from r ci cr approach is especially useful in recovering shapes of objects 
whose general structure is known but which may have little discernable texture in significant parts of their surfaces. A 
good example of such an object is the human face, where there is usually a significant amount of relatively untextured 
regions, especially when there is little facial hair, and where the general facial structure is known. Also described below 
is an application of the method to 3-D face modeling from muttiple images to generate a 3-D synthetic talking head. 

45 [0019] The talking head model comprises a collection of 3-D triangular facets, with nodes as vertices. The recovery 
is performed from a sequence of images taken with a camera having mostly unknown intrinsic and extrinsic parameters. 
[0020] In one embodiment a frontal image of the head, i.e., the face, is used as a reference image. Line-of-sight con- 
straints of 3-D facial nodes are imposed using the reference image. The 3-D model deformation is constrained by min- 
imizing an objective function that trades-off minimal change in local curvature and node position with fit to predicted 

so point correspondences and face appearance. 

[0021 ] For general camera motion with constant intrinsic parameters, three views are theoretically sufficient to recover 
structure, camera motion, and all five camera intrinsic parameters. For stability reasons, only one unknown intrinsic 
camera parameter is assumed, namely the focal length. The aspect ratio is assumed to be unity, the image skew to be 
insignificant, and the principal point to be coincident with the center of the image. 

System Overview 

[0022] Figure 1 shows an arrangement 100 including a camera 1 10, an image processing system 120, and a monitor 



3 



EP 0 907 144 A2 



130. The camera 110 can acquire a sequence of digital images 203-205 to be processed by the system 120. The 
sequence of images 203-205 can include various poses of an object for example, a head 10. Although two images 
should be sufficient to perform some type of 3-D recovery, *rt is suggested that a minimum of three images are used 
when the shape to be recovered is a face. The poses measured by the camera 1 10 produce a reference frontal image, 
s and left and right oblique views. An alternate arrangement is to have multiple cameras at predefined locations, each tak- 
ing a snapshot of the same face from a different angle. 

[0023] The system 120 can be a conventional computer system or workstation including input and output devices, for 
example, a mouse 121 and a keyboard 122. The system 120 includes one or more processors (P) 125. memories (M) 
126, and input/output (I/O) 127 interfaces connected by a bus 128. The memories 126 store data and instructions that 
io operate on the data. The data can be the digital images 2Ci ivc S. vie sequence, or "frames" acquired by the camera 
110. The system 120 can also be connected to a bulk storeys d&vics, for example, a disk 123. 

General Operation 

1 s [0024] During operation, the system 1 20 uses an appearance-based structure from motion technique to recover a 3- 
D model from initially unknown point correspondences in the images, and an approximate 3-D template. The model can 
be displayed on the monitor 130. The template can be a "wire-frame mesh templates" having a plurality of polygons, 
described in detail below. 

[0025] Figure 2 shows a flow diagram of a general method 200 for recovering a 3-D structure according to an embod- 
20 iment of the invention. The method 200 has an initialization phase 201 . and a refinement loop 202. Image registration 
210 is performed on images 203-205. Each image 204-205 is registered with a reference image 203. In the preferred 
embodiment, the image registration is spline-based as descr^J in U.S. Patent 5,611,00 issued to Szeliski on March 
11. 1997. 

[0026] In step 220, structure motion and parameters are io>jn>aisd using an iterative Levenberg-Marquardt batch 
25 approach. The appearance prediction is done using simple texture resampling. The predicted appearances 214-215 
corresponding to each of the images 204-205 are based on current image point correspondences and structure from 
motion estimates. The predicted appearances can then be used to refine the image registration during the loop 202. 
[0027] With the appearance-based structure from motion method according to the invention, initialization is first done 
by performing pair-wise spline-based image registration between the reference image and every other image in the 
so sequence. In other words, the sequence of images 203-205 acquired by the camera 1 1 0 includes at least two images. 
Better results can be obtained with a sequence of a greater number of images. 

[0028] The registration establishes a set of gross point correspondences for the image sequence. The camera param- 
eters and model shape are extracted from the point correspondences. Subsequently, the method iterates over the fol- 
lowing three major steps in the loop 202: 

35 

Appearance prediction 

[0029] In this step, for each image 204-205 other than the reference image, the appearance as constrained by the 3- 
D model 221. given the camera pose and intrinsic parameters, is determined and projected onto a hew image (214- 
40 21 5), i.e., predicted appearances. 

Spline-based image registration 

[0030] In step 21 1, each of the predicted appearances 214-215 is registered with the actual corresponding images 
45 204-205 to refine the point correspondences. Here, the loop 202 can terminate when a predetermined termination con- 
dition is reached, for example, a fixed number of iterations, or the average difference between the images 204-205 and 
the predicted appearances 214-21 5 is less than some threshold. 

Structure from motion 

so 

[0031] Using the refined point correspondences, estimate the new, and usually better estimates of the camera pose 
and intrinsic parameters, as well as the 3-D model shape. 

[0032] The use of appearance-based strategy is important because it accounts not only for occlusions, but also per- 
spective distortion due to changes in the pose of the object being imaged. In contrast to prior art techniques which use 
55 edges, or specific object features, e.g., the eyes, nose, and mouth, the entire predicted images are used here. 
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Tracking by spline-based registration 

[0033] In the preferred spline-based registration framework, each new image I2 (204-205) is registered to an initial 
referer.es image ^ (203) using a sum of squared differences formulation (1): 

5 

1 

10 " . - ... , . 

where the (Up vj)'s are estimates of optical flow on a per-pixel basis (x p y,). 

[0034] In Figure 3, spline control vertices {u p vj) are shown as circles (o) 310, and the pixel displacements {u it v>) 
are shown as pluses (+) 320. In the preferred registration technique, the flow estimates {Uf, vj) are represented using 
two-dimensional splines that are controlled by a smaller number of displacement estimates Qj and v $ which lie on a 
15 coarse spline control grid 330. This is in contrast to representing displacements as completely independent quantities 
having an underconstrained problem. The value for the displacement at a pixel / can be written as (2): 



or 



25 




30 where the Bj(x, y) are called the basis functions having non-zero over a small interval, i.e., there is finite support. Trie 
Wy = BjfXf, yj are called weights to emphasize that the displacements (u- p vj are known linear combinations of the 
control vertices (Qj, vj), 

[0035] In the preferred implementation, the spline control grid 330 is a regular subsampling of the pixel grid, e.g., 
x j =mx f , and y j = my t . Thus, each set of m x m pixels corresponds to a single spline patch 340. The bilinear basis 
35 functions for the spline can be expressed as: 

B/x, y) « max ((1- \x - Xjf/m) (1 -/y - y J / m), 0). 

[0036] Other bases are possible. The local spline-based optical flow parameters are recovered using a variant of the 
40 Levenberg-Marquardt iterative non-linear minimization technique. 

[0037] The spline-based equation (2) is modified to include the weights associated with a mask as follows (3): 

where m t j = / or 0 depending on whether the corresponding pixel is in the object or background area respectively. This 
is necessary to prevent registration of the background areas influencing registration of the projected model areas 
so across images. The value m f j can also be between 0 and 1 , especially during a hierarchical search where the images 
are subsampled and the intensities are averaged. 

General Structure from Motion 

55 [0038] The step 220 of recovering structure essentially involves trying to recover a set of 3-D structure parameters p,- 
and time-vi;! n'Ihq motion parameters 7} from a set of observed image features Uy. The general equation linking a 2-D 
image feature location u f j in frame / to its 3-D position p,-, is (4): 
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where /is the track index, and the perspective projection transformation PO is applied to a cascaded series of rigid 
transformation Tf^. Each transformation is in turn defined by (5): 



j 



where Pf® is a rotation matrix and tf$ is a translation applied after the rotation. Within each of the cascaded time-var- 
ying when the j subscript are present, and fixed when the j subscript is dropped. 
[0039] The general camera-centered perspective projection equation is (6): 



75 



where / is a product of the focal length of the camera and the pixel array scale factor, r is the image aspect ratio, a is 
the image skew, and (up, Vq) ' s *h e principal point, i.e., the point where the optical axis of the camera 110 intersects the 
20 image plane. In theory, for general camera motion with constant intrinsic parameters, three views are sufficient to 
recover structure, camera motion, and all five camera intrinsic parameters. For stability, only one intrinsic camera 
parameters matter is considered, namely the focal length, the aspect ratio is assumed to be unity. 
[0040] An alternative object-centered formulation can be expressed as (7): 

25 fx\ ( sx+r \*> . , "\ f sx ^ 
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rsy 



with a reasonable assumption that a - 0 and (u 0 , v 0 ) = (0, 0). Here, it is assumed that the (x, y z) coordinates before 
projection are with respect to a reference frame that has been displaced away from the camera by a distance f z along 
the optical axis, with s = f/t z and t| = 1A Z . It is possible to consider t 2 as the z component of the original global translation 
which is absorbed into the projection equation, and then set the third component of t to zero. 
35 [0041 ] The projection parameter s can be interpreted as a scale factor and r\ as a perspective distortion factor. The 
alternative perspective formulation (7) results in a more robust recovery of camera parameters under weak perspective, 
where r\ » 1, and assuming (uq, v 0 ) » (0, 0) and <j » 0, and P(x, y. zf » (sx, rsy) T . This is because s and r can be much 
more reliably recovered than tj, in comparison with formulation '£) where f and t z are highly correlated. 

40 Least-squares Minimization with Geometric Constraints 

[0042] The Levenberg-Marquardt algorithm is used to solve for the structure and motion parameters. However, with 
geometric constraints the method minimizes (8): 

where (9): 

I 1 

is the usual structure from motion objective function that minimizes deviation from observed point feature positions, and 
55 P() is given in (4) above, and (10): 



_ / T T m TxT 
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is the vector of structure and motion parameters which determine the image of point / in frame /. The vector a contains 
v ail of the unknown structure and motion parameters, including the 3-D points p,. the time-dependent motion parameters 

my, and the global motion/calibration parameters m g . The superscript T denotes a vector or matrix transpose. The 
weight c,y in (9) describes the confidence in measurement u,y. and is normally set to the inverse variance a,/ 2 . Here, the 
5 Cy can have a value proportional to the least amount of local texture indicated by the minimum eigenvalue of the local 
Hessian. 

[0043] The local Hessian H is given by: (1 1) 



10 



H = 



T I I T. I 2 
~w'x ' y^w' y_ 



where w is the local window centered at (x, y) and I y ) is the intensity gradient at pixel (x, y). If e mhi % is the minimum 
is eigenvalue at point / in frame y\ then: (12) 

c e min, ij 

9 max^-e^ 

20 . . . 

[0044] This is particularly important in the case of recovery of a face model because of the possible iack of texture on 
parts of the face, such as the cheeks and forehead areas. Using this metric for c,y minimizes the importance of points 
on these relatively untextured areas. To account for occlusions, eg is set to zero when the corresponding point is pre- 
dicted to be hidden. 
25 [0045] The other term in (8) is (13): 



^ eom ra;=X;ra/|A>r^/l 2 +P/!PrP/l 2 ^ 
/ 

30 



which is the additional geometric constraints that reduces the deformation of the template or reference 3-D model. The 
quantities with the superscript 0 refer to the reference 3-D model that is to be deformed, h t - is the perpendicular distance 
of point Pi to the plane passing through its nearest neighbors, here three. In other words, //,• is the best fit plane of the 
35 neighbor points of p /( and pxn , ~ d s is the equation of //,*, then (14): 

hi=Pfxn r d f 

a, is the weight associated to the preservation of local height, i.e.. to preserve curvature, and P/ is the weight associated 
40 with the preservation of the reference 3-D position. The weights can be made to vary from node to node, or made con- 
stant across all nodes. 

[0046] The Levenberg-Marquardt algorithm first forms the approximate Hessian matrix (15): 



where B($j) is a matrix which is zero everywhere except at the diagonal entries corresponding to the ith 3-D point. The 
so weighted gradient vector is (1 6): 

55 - 

where g, = (0 ... p', r ... 0) T f and (17): 
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from (14) and using the simplifying assumption that each node position is independent of its neighbors, although this is 
not strictly true. Here, e ff = u,j - f(a^) is the image plane error of point / in frame j. 

[0047] Given a current estimate of a, an increment 5a towards the local minimum is determined by solving (18): 

(A+\l)&a=-b t 

where X is a stabilizing factor which varies over time. A line-of-sight constraint is also imposed on the recovered 3-D 
point with respect to the reference image. 

Generating Predicted Appearance 

[0048] It is relatively easy to render the model given the facets and vertices of a 3-D surface model, described in detail 
below, e.g., the 3-D geometric constraints 221 of Figure 2. The object facets of the surface model are sorted in order of 
decreasing depth relative to the camera 1 10, and then rendered by texture-mapping the facets in the same decreasing 
depth order. The rendering technique can be any standard technique used in computer graphics. 
[0049] A 3-D model that is a good candidate for our proposed approach is the human face model. Its structure is 
known and using conventional stereo techniques are not very reliable because the human face usually has significant 
portions of relatively untextured regions. 

Mapping faces to a 3-D Computer Generated Face 

[0050] DECface is a system that facilitates the development cf applications requiring a real-time lip-synchronized syn- 
thetic talking head, see U.S. Patent 5,657,426, issued to Waters on August 12, 1997. 

[0051] DECface has been built with a simple interface protocol to support the development of face-related applica- 
tions. The fundamental components of DECface are software speech synthesis, audio-fae, and face modeling. 
[0052] Of particular importance to the present invention is the face modeling component It involves texture-mapping 
of frontal view face images, synthetic or real, onto a correctly-shaped wireframe. 

[0053] Topologies for facial synthesis are typically generated from explicit 3-D polygons. For simplicity, a simple 2-D 
representation of the full frontal view is constructed because, for the most part, personal interactions occur face-to-face. 
[0054] As shown in Figure 4, a wireframe model 400 is superimposed on a facial image. The model 400 includes 200 
polygons of which 50 represent the mouth and an additional 20 represent the teeth. The jaw nodes are moved vertically 
as a function of displacement of the corners of the mouth. The lower teeth are displaced along with the lower jaw. Eye- 
lids are created from a double set of nodes describing theupper lid, such that as they move, the lids close. 
[0055] The canonical representation is originally mapped onto the individual's image mostly by hand. This requires 
the careful placement of key nodes to certain locations, the scvners of the lips and eyes, the placement of the chin and 
eyebrows, as well as the overall margins of the face. 

Mapping Faces using one input image 

[0056] As mentioned in the previous section, mapping new faces to DECface involves texture-mapping frontal view 
face images (synthetic or real) onto a correctly-shaped wireframe. The original method to generate DECface with a new 
face is to manually adjust every node, which is a very tedious process. A "generic" separate face whose DECface topol- 
ogy and 3-D distribution is known is used as a reference during the process of moving each node within the new face 
image. This node-moving process is equivalent to the transfer of z information from the generic face to the new face. 
Methods to automate this process by using templates of facial features such as the eyes, mouth, and face profile have 
also been used For a detailed description please see, EP application 98304034.6, "Automated Mapping of Facial 
Image to Wireframe Topology", filed by Digital Equipment Corporation, inventors being Kang and Waters. 
[0057] Because only one face input image is used, to gens:£!3 the appropriate 3-D version of DECface, the canonical 
height distribution is preserved. This is, however, not always desirable, especially since many human faces have signif- 
icantly different facial shapes. As a result, to preserve as much as possible the correct shape, we use three input 
images, each showing a different pose of the face, with one showing the frontal face pose. It is posstole, of course, to 
use more than three images to achieve the same goal. 
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Mapping faces using three input images 

[0058] In the preferred embodiment, three images of a face at different orientations are used. An image of a frontal 
pose is used as the reference image. As before, the camera parameters, intrinsic and extrinsic, are generally not known. 
5 For simplicity, one can assume that the aspect ratio is one, the image skew is zero, and the principal point is at the 
image center. The point correspondences between the generic face and the reference face has been done as in 
described above. This is the same as assuming that the reference shape of the model has been initialized. Note, how- 
ever, that the point correspondences across the image sequence are not known. 

[0059] The values a,- and 0, in (13) are set to 0.25. As mentioned above, the feature track fine-tuning step involves 
io us. spline-based tracker on the predicted appearance and actual image. However, because the prediction does 

not involve the background, only the predicted face image portion of the image is involved; the weights associated with 
the background are set to zero in the spline-based tracker. 

[0060] The method may fail when the change in the appearance of the object is too drastic from one frame to another 
in a sequence. In an 3-D face modeling application, rotations of up to about 15° between images are well tolerated. 
15 [0061 ] A variant of the method would involve the direct incorporation of the optic flow term into the objective function 
(8) to yield (19): 

e a //^= E srm( a ) + * geom( a )+* flow( u ) 

20 where (20): 

25 

with lj(U{j) being the intensity (or color) at Uj on frame /, and y f y is the weight associated with the point u%. Note, that in 
this particular application of facial model recovery, the value of un is kept constant throughout because the first frame 
is the reference frame. 

30 [0062] One problem with directly embedding this term in the structure from motion module is that the flow error term 
is local and thus unable to account for large motions. It would either require that the initial model pose be quite close to 
the true model pose, or the addition of a hierarchical scheme similar to that implemented in the spline-based registration 
method. Otherwise, the method is likely to have better convergence properties when the tracking is performed outside 
the structure from motion loop. In the present implementation, while having the small perturbations of the model pose 

35 would be desirable from the computational point of view, although not from the accuracy point of view, this is not a 
requirement. 

[0063] In addition, using the flow error term directly may not be efficient from the computational point of view. This is 
because at every iteration and incremental step, a new predicted appearance has to be determined. This operation is 
rather computationally expensive, especially when the size of the projected model is large. Having the tracking module 
40 only loosely coupled with structure from motion results in fewer number of iterations in computing the predicted object 
appearance. Finally, there is the non-trivial question of assigning the weights jy relative to the structure from motion and 
geometric constraint related weights. 

[0064] Geometric constraints on the face deformation in other forms can also be used. An example would be to use 
the most dominant few deformation vectors based on Singular Value Decomposition (SVD) analysis of multiple training 
45 3-D faces. A similar approach would be to apply nodal analysis on the multiple training 3-D faces to extract common 
and permissible deformations in terms of nonrigid modes. 

Summary 

. so [0065] Described is an appearance-based structure from motion method that enables the direct extraction of 3-D 
models from a sequence of uncalibrated images. It is not necessary to precompute feature correspondences across the 
image sequence. The method dynamically determines the feature correspondences, estimates the structure and cam- 
era motion, and uses this information to pre^'-c: the object appearance in order to refine the feature correspondences. 
[0066] It is understood that the above-descrifcsti embodiments are simply illustrative of the principles of the invention. 
55 Various other modifications and changes may be made by those skilled in the art. which will embody the principles of 
the invention and fall within the scope thereof. 
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Claims 

1. A computerized method for extracting a three-dimensional configuration information from a sequence of two- 
dimensional images of an object, comprising the steps of: 

5 

(a) registering each image in the sequence with a reference image to determine image features; 

(b) recovering structure and motion parameters from the image features using geometric constraints; 

(c) generating a predicted appearance for each image using the recovered structure and motion parameters; 

(d) registering each predicted appearance with the corresponding image; and 

10 repeating the recovering (b), generating (c) and registering (d) steps until a termination condition is reached. 

2. The method of claim 1 . wherein the registering is done using spline-based image registration. 

3. The method of daim 1 , wherein the geometric constraints are imposed directly on the 3-D wireframe model of the 
15 object. 

4. The method of claim 1 , wherein the object is a face. 

5. The method of claim 1 , wherein the sequence of images includes at least two images. 

20 

6. The method of claim 1, wherein the repeating terminates when an average pixel value difference between each 
predicted appearance and each corresponding image is less than a predetermined threshold. 

7. The method of claim 1 , wherein camera parameters are generally unknown. 

25 

8. The method of claim 3, wherein the wireframe model is composed of triangular facets. 
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